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(N| 1 I. INTRODUCTION 

> 

In order to obtain higher accuracy in the templates for analyzing gravitational waves from binary black holes, 
gravitational spin-interaction terms have to be taken into account beyond the leading order ones which are of the 
formal order 1/c 2 (in this counting, spins are treated of the order zero in terms of 1/c 2 ), where c denotes the speed of 
light. At the formal order of 1/c 4 , several Hamiltonians or Lagrangians have been determined already: The spin-orbit 
coupling ones, Hgo = Hs P + Hg p 3, are given in also see [2j, the spin(l)-spin(2) coupling ones, Hs 1 s 2 + Hg 1 g 2p 2, can 
be found in U, also see [|], and the following Hamiltonians H S 2 Sipi +H S 3 pi +H S 3 p2 +H S 2g 2P2 , H S 2 S 2, H SlS s +H S2S 3 
have been derived in recently, the dynamics corresponding to the Hamiltonians H S 2 +H S 2 +H S 2 p2 +H S 2 p 2 has been 
announced, In this paper the missing S^p-Hamiltonians Hsfp ± + Hs'^p 2 and Hg2 S2Pi +Hg2 Sipi> are calculated in the 
(3jT)' Arnowitt-Deser-Misner (ADM) canonical formalism by applying the space-asymptotic Poincare algebra, 0,0- These 
i— i , expressions contain post-Newtonian leading order quadrupole deformation effects represented by spin-squared terms 
(the coefficients in the quadrupole terms tell that the treated bodies are black holes, see, e.g. @). Because of the 
power of the Poincare algebra in controlling the higher post-Newtonian dynamics we are also able to give reasonable 
qq ■ expressions for the spin-squared S ,2 p 2 -Hamiltonians to linear order in G, i.e. H S 2 p2 +H S 2 p2 . Unfortunately, only after 
completion of the calculation of all static G 2 -terms Hg2 a comparison with the result presented in can be made 
CN ■ because of quite different approaches. 

Our units are very often c = 1 and also 167rG = 1, where G is the Newtonian gravitational constant. Greek indices 
Q\ • will run over 0, 1, 2, 3, Latin from the beginning of the alphabet over 1, 2 (denoting black-hole numbers) and from its 
midst over 1,2,3. For the signature of spacetime we choose +2. r = r a \ > = \x\ — x\\ will denote the Euclidean distance 
OO , between black hole a and b, and r n l = r ab n l ab — x\ — x\, {n l ab ) = n ab = n. 

> '■ 

II. THE POINCARE INVARIANCE 

The starting point for our approach will be the Poincare algebra. The Hamiltonians we wish to calculate have to 
fulfill this algebra using standard Poisson brackets for the fundamenal position, linear momentum and spin variables, 
x aiPai, (Pai) = P, and S a i, (S a i) = S, respectively, so the plan is to make the most general ansatz for these Hamil- 
tonians and plug them into the algebra to see how restrictive it is upon them. Of course, the Hamiltonians cannot 
be determined uniquely because all Hamiltonians which are canonically equivalent fulfill the algebra, so there will be 
degrees of freedom left which can be fixed only by choosing an appropriate representation. We will choose to work 
within the ADM canonical formalism using generalized isotropic coordinates, Q. With the aid of a reasonable ansatz 
for the source terms in the constraint equations we are able to fix all the coefficients of the Hamiltonians in question. 

The Poisson bracket relations defining the Poincare algebra read, see, e.g. @, 

{P i ,P j } = 0, {P 4 ,ff} = 0, {J h H} = 0, (2.1) 



{Ji, Pj} — eijkPk, {Ji, Jj} — tijkJk, (2.2) 
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{Ji, Gj} = f-ijkGk, (2.3) 
{Gi,H} = Pi, (2.4) 

{Gi, Pj} = -^H5 i: j, (2.5) 

{Gi,Gj} — —-^tijkJk- (2.6) 

The total linear momentum Pj and the total angular momentum Ji are respectively given by Pj = 5Z a p a i an< ^ 
Jj = ^oi) where J a i = eijkX J a Pak + S a i- The crucial equation which needs to be checked and which finally will be 
used for determining the new Hamiltonians is Eq. (|2.4|) . It describes the time evolution of the center-of-mass vector 
Gi. In the post-Newtonian context, the importance of Eq. (|2.4[) for the fixation of the Hamiltonian has been stressed 
already in Ref. [l(|, which is the source paper for Q. Obviously, the Hamiltonian in Eq. (2.5) needs to be known 
to one order in 1/c 2 less compared with Eq. (|2.4[) . and Eq. (2.1) only tells that the Hamiltonian has to be invariant 
against 3-dimensional translations and rotations. The second post-Newtonian Hamiltonian (2PN-Hamiltonian) of 
formal order 1/c 4 which enters the equations above reads 

H = Hn + Hipn + H 2 pn + HgQ N + HgQ N + H S 2 + H s s p + H S 2 p 2 + H S i (2.7) 

with 

H S i = H Sl s 2 + H S 2 + H S 2 . (2.8) 

The Hamiltonians Hn, Hipn and H2PN are the point particle ones given in e.g. @. The Hamiltonians H$q N , Hs 1 s 2 
and H S 2 + H S 2 , as far as being of order 1PN, were (re-) calculated in [f|, also see 

[llj], E3- Th e Hamiltonian H?"'~ 
is known from |l[ and HsiS2(p 2 +G) has been derived in [3[. The Hamiltonian Hg3 p splits into the parts 



2PN 
SO 



H S 3 P - H s s pi + H S 3 p2 + H S 3 pi + H S 3 p2 + H SlS 2 pi + H S 2 S2pi + H S2S 2 p2 + H S 2 Sip2 . (2.9) 
and Hg2 p 2 is an abbreviation for 

H S 2 p 2 = H SlS2p 2 + H S 2 p 2 + Hg2 p 2. ( 2 -10) 

The contributions Hg3 pi , H S 3 p2 , Hg 2 g2 p2 and H SlS 2 pi were all calculated in [5J] (Sb = m &<!;,), 

H slp2 = -GS v (n u x p 2 ) (4- - , (2 .ii) 

\ r 12 r 12 / 

H s s pi = H s s p2 (1^2), (2.12) 
„ „3m 2 f 3alp 2 -(S 2 x n 12 ) 6(oi«ni 2 ) p 2 - (S 2 x ax) 15(a 1 -ni 2 ) 2 p 2 -(S 2 x n 12 )\ 

H S 2 S 2 lP 2 = ~ G 7Z- I 14 + ZA 14 ) > \ 2 - 13 > 



Ami 



H Sl si Pl = H S2 sf P2 (1^2). (2.14) 



3 



The Hamiltonian Hgt splits into the parts 



Hs l — H S 2 S 2 + H SlS 3 + H S 3g 2 + H S i + H s i 



(2.15) 



The Hamiltonians H S 2 S 2 + H SiS z + H S 3 g 2 are given in Q and the Hamiltonians H S i + H S 4 will be shown to be zero 
in Sec. VI below. The remaining Hamiltonians will be chosen in the form, with to be determined /i-coefficients, 



G 



Si ■ (ni2 x p x ) (jiiSl + [l 2 (Si ■ n 12 ) 



(2.16) 



H sip2 = H sfpi (1 <- 2) , 



(2.17) 



H s 2 1 s 2Pl -~r 

'12 



\xzS\Si ■ (ni2 x p x ) + fa (Si ■ nu) S 2 ■ (Si x p x ) + // 5 (Si ■ m 2 ) S 2 ■ (n i2 x p x ) 
+ (Si • S 2 ) Si • (n i2 x p x ) + ^ 7 (5i • ni 2 ) (S 2 ■ n X2 ) Si • (n 12 x p t ) 
+ f»i 2 ■ (Si x S 2 ) (a*sSi ■ p x + (j, 9 (Si ■ nu) (Pi ■ n 12 ) 



(2.18) 



#s»s lPa = H sjS2Pl (1 <- 2) . 



(2.19) 



Notice that not all terms are independent of each other. They are connected via the following identity for the mixed 
product of three vectors in three dimensions, 



(Ui ,U 2 ,U 3 )U = (UUi) U 2 x U 3 + (UU 2 ) U 3 xUi + (UU 3 ) U x x U 2 
= (U,U 2 , U 3 ) Ui + (Ui,U, U 3 ) U 2 + (Ui ,U 2 ,U) Us 



(2.20) 
(2.21) 



with definitions (UUi) =U-U\ and (Ui , U 2 , Us) = Ui ■ (U 2 x Us)- This will help us later to understand that some of 
the /J coefficients must be zero, because some terms can be shifted into others. The ansatz for the S 2 p 2 Hamiltonian 
reads, using a, /3, and 7 coefficients, 



'12 



Qi(pi-Si) + a 2 p\S\ + as (Pi ■ ni 2 ) S\ + a±p\ (Si ■ n i2 ) 
+ a 5 (p 1 ■ ni 2 f (Si ■ ni 2 f + a 6 (p 1 ■ n 12 ) (Si ■ n i2 ) (p x • Si) 



Pi (Pa ' Si) 2 + P 2P 2 2 S 2 i + (3s (P 2 • ni 2 f S\ + /3 4 p| (S x ■ n 12 ) 
■ 05 (Pa • "12) 2 (Si • ni 2 f + 6 (p 2 • m 2 ) (Si • n i2 ) (p 2 • Si) 



(2.22) 



-7i (Pi -P 2 )Sl+j 2 (pi -p 2 ) (Si ■ ni 2 f +73 (p x • Si) (p 2 • Si) 
- 74 (Pi • n i2 ) (p 2 • Si) (Si • n i2 ) + 75 (p 2 • nia) (pi • Si) (Si ■ n i2 ) 



- 76 (Pi • "12) (p 2 ' "12) Si + 77 (p! ■ ni 2 ) (p 2 • n i2 ) (5i • ni 2 )' 



i? S 2 p2 (1 «-» 2) 



(2.23) 



The center-of-mass vector G which enters in Eq. (12. 4|) is given by 
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G — Gn + Gipn + G2PN + G 1 ^) + Glo^ + Gs 1 s 2 + Gsf + ^SJ- (2.24) 

The point particle expressions Gn, Gipn and G 2 pn are given in The parts Gg™ and G 2 S ^ N are presented in 
[l[ and GjjS 2 was calculated in 0]. What is left is to give the expression for G S 2 and G52. Here in this section we 
will also make a general gauge invariant ansatz for them to show how they are involved in the algebra. It reads with 
coefficients v\, v%, V3, i>a and v§ 



G S 2 = G 



m 2 
mi 



(Si.nia)Si , (Sx-n 12 ) 



V3X2) 



(2.25) 



plus the expression with (1 <-> 2). 

Later we will calculate G S 2 directly by an appropriate covariant source for the Hamilton Constraint, which will 
fixate all the ^-coefficients. 

III. FULFILLMENT OF THE POINCARE ALGEBRA 

We concentrate on Eq. (|2.4p which gives rise to the following equations for the coefficients: 



Ml 

t'2 
M3 
M4 
M5 
M6 
M7 



m 2 
4m? ' 

5 ma 
~4m?' 



15 
2m? 
15 



Ms, 



+ M9 +Ms, 



M9j 



2m? 
15 



Ms, 



2m? 



M9, 



(3.1) 
(3.2) 
(3.3) 
(3.4) 
(3.5) 
(3.6) 
(3.7) 
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3 m 2 m 2 m 2 
" 2 = ^ - 2^71 + ^3-5, (3.9) 



9m 2 m 2 3m 2 ,„ in , 

"3 = -q-3 - 77— 76- o^3^, (3.10) 



8 2mi 2 m? 

9 m 2 m 2 m 2 
i ra\ 2ra\ ^ 2 2m\ 



» ni 2 ni 2 111,2 . . 

"4 = -o33-— 72-7^3^2, (3.11) 



15 m 2 m 2 5m 2 
"5 = --T-z -^T-77 + ^-3^2, (3.12) 



4 to? 2mi ' 2 m? 
15 m 2 m 2 3m 2 

-j ? 74 H q- 

4 mi mj 



io 7/12 "'2 . 0"*2 ,„ 1 „n 

a 6 = 77 74 H o-vi, (3.13) 



A = ^-f^73 + ^, (3.14) 
wi\m 2 2m 2 2m\m 2 

P2 = -z 71 - o ' 3 - 15 ) 

m\m 2 2m 2 2m\m 2 

^ 3 4mim 2 2m 2 ^ 6 2m,\m 2 5 ' ^ ^ 

P4 = -: » 72 + ^ , (3.17) 

4mim 2 2m 2 2m,\m 2 

= —yr-^^—^, (3.18) 
a 3 mi 2^ 2 

Pe = 74 + , (3.19) 

mim 2 m 2 mim 2 

75 = 74 2 j, (3.20) 



mj mf 



3 



^3 = 2-^2, (3.21) 

^4 = (3-22) 

The Poincare algebra obviously fixes the Hamiltonian H S 3 pi (and thus also H S 3 p2 ). It also restricts five from seven 
coefficients for H S 2 S2pi (and thus also for H S 2 Sip2 ), and two out of five coefficients for the center-of-mass vector G S 2. 
At this stage of our investigations, the coefficients in H S 2 p2 and G52 can still be changed by the application of a 
canonical transformation generated by 

G ( 

9slp =— s l (°"iPi ' n i2 + ^2^2 • n i2) + Si ■ n 12 (03S1 • p 1 + a 4 Si ■ p 2 ) 

r !2 V , 
V (3.23) 

+ (Si ■ n l2 f (a 5Pl ■ n l2 + a 6 p 2 ■ n 12 ) J 

with coefficients c^, fc = 1, 2, 3, 4, 5, 6. Later on we shall see that two coefficients of the center-of-mass vector G S 2 , and 
thus many other coefficients which are connected with them, can uniquely only be fixed with the aid of the explicit 
expression for the energy density of a static source. Coming back to the Poincare algebra, all its other equations are 
trivially fulfilled and impose no further restrictions on our coefficients. 



IV. THE ADM GENERALIZED ISOTROPIC COORDINATES REPRESENTATION 



To get a unique representation of the Hamiltonians, we now have to fix the coordinates. We will use the ADM 
formalism and generalized isotropic coordinates. This makes it very easy to calculate interaction Hamiltonians if 
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the source expressions in the constraint equations are known. Before imposing constraint equations and coordinate 
conditions, the Hamiltonian in the ADM framework reads [?J [H, 



H = / d 3 x{NH - N*7U) + E [ 7ij ] 



(4.1) 



where respectively N and N l denote lapse and shift function, which are merely Lagrange multipliers. The super- 
Hamiltonian Ti and the supermomentum TLi densities decompose into gravitational field and matter parts as follows: 



H = H 



field 



field 



where the field parts are given by 



^field ^_ 



7R+2 hijK ij ) 2 



7~U = 7~C 



(4.2) 



(4.3) 



Here 7 is the determinant of the 3- metric 7 y = gij of the spacelike hypersurfaces t = const., whereas the determinant 
of the 4-dimensional metric g^ v will be denoted g. The canonical conjugate to 7^ is 7r y . R is the Ricci scalar of 
the spacelike hypersurfaces and ; denotes the 3-dimensional covariant derivative. The expressions for lapse and shift 
functions in terms of the metric are N = {— g 00 ) -1 / 2 and N l — j^goj, where 7 U is the inverse 3-metric. As in Q, we 
will assume that the relation between field momentum -k ij and extrinsic curvature Ka is the same as in the vacuum 



tt« =-V7(7 i V-7 tf 7 W )^w- 
The energy E is defined by (herein, the comma "," denotes partial 3-dimensional derivative) 



E = d> d 2 s,( 7 j 



(4.4) 



(4.5) 



E turns into the ADM Hamiltonian Hadm after imposing the constraint equations H = Hi — and appropriate 
coordinate conditions. Comparing the constraint equations with the Einstein field equations, projected onto the 
spacelike hypersurfaces, results in 



(4.6) 



where ^J—gT^ 11 is the stress-energy tensor density of the matter system. The timelike unit 4-vector n M = (— JV, 0, 0, 0) 
points orthogonal to the spacelike hypersurfaces. 

The generalized isotropic coordinates 0] , also called the ADMTT gauge, are the most often used and best adapted 
coordinate conditions for explicit calculations and they are defined by 



Hi = + h]? 



with ip 



1 



or 37^ 



7mV 



, and 7T l 



0, 



TT 13 + 7T 



TT ' 



(4.7) 



hlj has the properties ftM = — (transverse and traceless) and the same holds for 7Trf T which is the canonically 
conjugate to hJ T . ff lJ denotes the longitudinal part of 



The ADM Hamiltonian and the center-of-mass vector read 

E = H ADM = - [ d 3 xAcj), 



G, 



<Fxx l A<j), 



(4.8) 



where (j> is expressed in terms of the canonical matter variables and the canonical field variables hj? and 7r^ T of the 
field-reduced phase space. Herewith the Poincare algebra can be verified. The most elegant way for the obtention 
of an autonomous conservative Hamiltonian in the matter variables is via the Routhian approach, see, e.g. [l3l |. 
However, to the order of our investigations 7r^ T does not play any role so we may stay on the Hamiltonian level for 
the field variables. 

The constraint equations explicitly read 
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n " ;fc 



/T; 



(4.9) 
(4.10) 



To include momentum squared terms interacting with spin-squared terms we have to make an ansatz for the sources 
generalizing the point particle terms without spin. The most general ansatz for our purpose reads for the source in 
the Hamilton constraint, cf. 



E 



- ^0%^ - \Pb ■ K x 9) + (j ij Pbi p bj + ml) 1 '* + X^Qfdid. 



1/2 



Pt 



— (Pb ■ 9)Q l b J Pbidj + —al(p b ■ df - \ 8 p b ■ (a b x d) Q^d^ 



(4.11) 



with the quadrupole tensor 



Ql 3 = « - 3«'% . (4-12) 

where by definition a b = S b /m b holds and 6b = 6{x l — x\) with J d 3 x 6 b = 1 and d = {di),di = A. Notice that 
purely Laplacian source terms, A6b, apart from the ones appearingvia Q 13 , are left out because they would lead to a 
distributional Hamiltonian (cf. paragraph on Breit's equation in [14j |) which is of no interest here. 

The ansatz for the source of the momentum constraint reads, also not taking into account terms of the form A6 b 
apart from the ones appearing with Q lJ (taking into account those terms would lead to no effects on our Hamiltonians 
when controlled by Poincare algebra), cf. [5J], 



Qb l ( XbPbkdidi + \ 6 pbidkdi + \ 7 (p b ■ 9) 8ud k I + \ia 2 b (p b ■ d)d. 



1 



(a b x9),(l- -Q k b l d k d t ) - -pb 



1 



(4.13) 



This ansatz is not 3-dimensional general covariant, but it is of sufficient general form that will lead to all the 
searched for contributions of the Hg3 p and Hg2 p 2 Hamiltonians. It also correctly holds J d 3 x 7Y™ atter = P t and 
/ d 3 x ey-fc£ J ?i™ atter = Ji. To the order needed for a 2PN-Hamiltonian for self-spin interaction, the Hamilton con- 
straint expands as 



-A</> (2) = UT 2 



matter 

(2) : 



A J. T/mattcr 

-A0( 6) = n { 



He) 



-a</» (8) = n 



matter 
(8) 



-A0 (4 ) = H 



(4) 



-H 



(2) 



(4) 



/oyrnatter ± 
8^(6) <?(2) 



n m 0( 4) ) + g^«(2) 



H 



(2) 

matter 



matter ±2 



2<P(2), l J ft (4) lJ 



/matter j 3 



512 "(2) ^(2) 



(4) V(4) T "-(2) 



\ I /a/ matter 1 2 

(6)) 64 ( W ^(2) 



(2) 



121 1 •«,;„ 1 " *'"(3)'"(5) ~ ^(2)^(2)^^ + -{hj^ ijk ) 2 



'{2)9(4)1 
(td), 



(4.14) 
(4.15) 

(4.16) 



( "td" means total derivative) and the momentum constraint can be written as 



~ 13 -T/mattcr 

Z-V 1-ivmatter l± zM \ 



(4.17) 
(4.18) 



The integral over — A0( 8 ) is the one leading to the new Hamiltonians to linear order in G. The hj^ term will drop 
out from the calculation of the Hamiltonians in question. The expressions needed for this integral read 
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a/mattcr 
«(2) 


= ^m b <5h, 

6 






a/mattcr 


= E 

6 






5b, 


o/matter 
"(6) 


= E 

6 




A 8 p b • (a 6 x 9) Ql J didj - 





a/matter 
"(8) 



(3), 



a/mattcr 
«(5)i 



m 6 



+ — (P b ■ d)Ql'p bi dj + —a 2 b ( Pb ■ df ) S b , 



^ \ V {i 'm b 64^ (2) m 6 

b v 

^ mi - -m 6 (at x 9). J £ t 



Qb' ( ^bPbkdidi + X 6 pbid k di + A 7 (p 6 • d) 8 u d h 



(4.19) 
(4.20) 



(4.21) 
(4.22) 
(4.23) 
(4.24) 
(4.25) 



and 



'(6) 



-II 

6 (6) 

D (6) 
& (6) 



n 



4>w = 4G E 



m b 
~2~ 



2m 6 2 ?"b „ ^ab^f) 



a^b 



A — laymattcr 



(4.26) 
(4.27) 



4G V ( A 8Pb • (a b x 9 b ) Qf£> w % + \i^-Q l b d bi d bj + ^( Pb ■ d b )Q^ Pbl d b] + —a 2 b (p b ■ d b ) 2 }- 
^ \ <tm b m b m b I r b 



a#b 



-A" 



m b r ab r b 



" t T_/mattcr ± 

-r(« (4) 0(2) 



(2) 



>(4) 



^ I M H (2) ^(2) 



-A' 



(4.28) 

(4.29) 
(4.30) 
(4.31) 



-27 



©1 J 



o_/matt< 
2«(3)fe 



(5)fc 



(0(2)7f*3)).i 



with 



(4.32) 
(4.33) 
(4.34) 



We do not need to know fr^ and ff^ in full detail as we can apply partial integrations within our used analytical 
regularization procedure. The integrals are calculated with the aid of the same methods and rcgularization formulas 
as outlined in, e.g. The results are (the symbol ~ indicates that only the relevant contributions are given) 
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H m = [n^d?x^GY J \ r ^PlQ l id ai d a ~, (4.35) 

J _lu b < ab 

#(8)2 = J —H^^x 

^ ~ ^2 ma ( X »Pb ' ( a b x 9 b) Q l b J d bi d bj + \i-^-Q b d bi d bj + — (p b ■ d b )Q l b J p bi d bj 

+ ^ b al{Pb - 9b)2 )h (4 - 36) 

^8)3 = / -\n?^^ x = / -\n^ tter (^6) + <A(6) + C) + C)) d3x > ( 4 - 3? ) 

#(8)3 = #(8)2, (4-38) 

HH )3 = 0(G 2 ), (4.39) 

H% 3 = 0(G 2 ), (4.40) 

#(8)3 = Jl^(2)^ ) d 3 x = |-^ (2) (^ 3) ) 2 d 3 x, (4.41) 

#(8)4 = |-^f 4 f tor 0( 4 )d 3 x (4.42) 

- i G E ^U^- + O (G 2 ) , (4.43) 

a^b 

#(8)5 * / -^#^ tter ^ 2 )d 3 ^ - O (G 2 ) , (4.44) 

#(8)8 = / ^(2) (^ 3) ) 2 d 3 x = - JTg 3) (4.45) 
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# (8)9 = / 27r^ 5) d 3 z 



5-3- f A 4 - -A 5 - -A 7 



mm 2 
GA 5 



6S 2 (Pi • P 2 ) - !8S 2 (pi • ni 2 ) (p 2 ' "12) 



+ 



+ 



3 3 

6 (Si • Pl ) (Si • ft) - 18 (Si • ni 2 ) (p 2 ■ n 12 ) (Si • Pl ) 



15 ( Pl ■ n 12 ) {p 2 ■ n 12 ) (Si • ni 2 ) 2 - 21 ( Pl ■ p 2 ) (Si • n 12 f 



+ 



m\r\ 2 

GX e 
m\r\ 2 

+ 6 (p 2 ■ n 12 ) (Si • n 12 ) (p 1 • Si) + 6 (p 1 ■ n 12 ) (Si • n 12 ) (p 2 ■ Si 
- 2 ( Pl ■ Si) (p 2 • Si) + 3 ( Pl • m 2 ) (p 2 • m 2 ) S 2 + 7 ( Pl • pa) S 2 
GA 7 



mfrJa 



15 ( Pl ■ n 12 ) (p 2 • m 2 ) (Si • n i2 ) 2 + 3 ( Pl • p 2 ) (Si • n i2 ) 2 



+ 6 (p 2 • m 2 ) (Si • n 12 ) (p 1 • Si) - 18 ( Pl • m 2 ) (Si • m 2 ) (p 2 • Si) 
+ 6 ( Pl • Si) (p 2 • Si) + 3 ( Pl • m 2 ) (p 2 • ni 2 ) S 2 - ( Pl • p 2 ) S 2 
GA 6 



+ 



+ 



3m 2 rf 2 

GA 6 / 
m 2 rf 2 I 

m\v\ 2 \ 



(4.46) 



- 2S X ( Pl • p 2 ) + 6S-L ( Pl • n i2 ) (p 2 • n i2 ) 



- 60 (Si • n i2 ) 2 1112 ■ (Pi x S 2 ) + 12S 2 ni 2 • ( Pl x S 2 ) + 24 (Si • m 2 ) Pl • (S 2 x Si) 



60 ( Pl • m 2 ) (Si • n 12 ) n 12 • (Si x S 2 ) + 12 ( Pl • Si) n 12 • (Si x S 2 ) 



+ 12 (Si -ni 2 ) Pl - (Si x S 2 ) 



+ (l<-»2). 
We conclude 



ff , rr 3 G 2 / 3(Sl-m 2 ) 2 -Sr , 

#( 8) i + #(8)4 = ~ A —P 2 I ] + (1 - 2) , 



#(8)3 + #(8)8 - 0. 



' 12 



#(8)2 + #(8)3 - ~G 



m 2 



Ai 



P 2 3(Si • m 2 ) 2 - S 2 



+ ( A 3 - ^A 2 ) St 



3(Pi -n 12 ) 2 -pi 



+A 2 



-A 8 



3( Pl • n i2 )(Si • ni 2 )( Pl • Si) - ( Pl • Si) 2 
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(-3S 2 + i5(Si-m 2 ) 2 ) Pl .(Si xm 2 ) 



r\ 2 



+ (!<-> 2), 



The coefficient equations following from our considerations read 



(4.47) 



(4.48) 
(4.49) 
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m 2 m 2 f, 2 \ n m 2 1 



ai = rr^ 2 ' " 2 = 2^f + 2A 3 - 3 A 2j > "3 = -3^ (^A 3 - -A 

3 m2 m,2 
q 4 = — -Ai — o, «5 = 0, a 6 = -3A 2 — 2' 
2 mf mj 

A = ft - ft = ft = 0, ft: 



\2 



4mim2 ' 4mim2 ' 

1 / 23 

71 = — o 6A4 - 2A 5 - 3A 7 + — A 6 

\ 3 

72 = — (-21A 6 + 3A 7 ) , 

73 = — (6A5 - 2A 6 + 6A 7 ) , 

74 = —9 (6A 6 - I8A7) , 

777, J 

75 = -4 (-I8A5 + 6A 6 + 6A7) , 

76 = — (-I8A4 + 6A 5 + 9A 7 + A 6 ) , 

77 = — a- (-15A 6 - 15A 7 ) , 

mf 

m 2 m 2 12 

Ml — — 3A8, /i2 — 10 — 3^8, M3 — — jAq, 

mf m i m i 

1 , -60 12 60 

M4 — — 2 (24A 6 — 12A 7 ) , /i 5 — 5-A6, M8 — — 2^7, M9 — — 2~^7, 

mf mf mf mf 

Me = M7 = 0. 



Now we make use of the Eqs. (|3.1[) to fl3.22[) and end up with the solution for all the coefficients still having 2 degrees 
of freedom left parametrized by v\ and 



m 2 ( 7 A m 2 / 5 A m 2 / 27 A 3 m 2 

mf \ 4 / mf \4 / mf \ 8 / 8 mf 

m 2 f 21 



"6 = — 3 \ —r + 3^i 

mf V 4 



ft=ft=ft=ft = 0, ft = -- , (3, 



4mim 2 ' ^ 4 4mim 2 ' 
I 9 1 ,„ 3 1 



7i = T^-fC- 1 -^), 72 = -7-2, 73 = 7-3(2 + ^1), 74 = -7^—2: 75 = — 2 (-3 - 3n) 



1/9 „ A 15 



m2 5 m 2 3 3 15 3 15 

Ml = 4mf M2 = ~4r< M3 = 2mf ' ^ 4 = 2mf ' M5 = ~2mf M8 = 2r< ^ 9 = _ 2mf ' 
Me = M7 = , 

3 3 1 

V2 = -, ^3 = J, ^4 = -5 - ^5 



Ai - -, A 2 = - i -^ 1) \ 3 = --- + v 5 , ^ = "^+18-6-. As -I + y 
A7 = g, A 8 = --. 



4.50) 

4.51) 

4.52) 

4.53) 

4.54) 

4.55) 

4.56) 

4.57) 

4.58) 

4.59) 

4.60) 

4.61) 
4.62) 



4.63) 

4.64) 

4.65) 

4.66) 

4.67) 

4.68) 
4.69) 
4.70) 

4.71) 
4.72) 



12 



It is worthy to point out, that the free parameters v\ and z/5 are related to a canonical transformation which is 
given by the generator 

Ss?pi = —Pr {-"5 S i(Pi ■ "12) + ^(^i ■ "i2)(Pi • Si)) . (4.73) 
m 1 r 12 

Refering to Eq. (|3.23p this means that we may choose u\ = —^^5, 03 = ^%v\, and er 2 — 04 = 05 = o§ = 

to generate the terms related with these coefficients. So all Hamiltonians parametrized by v\ and 1/5 are canonically 
equivalent and we are free to give them any value. To fixate them we need to calculate Gg 2 explicitly. We use Eq. 
(14. 8|) and as a matter source for the static case we adopt the result from [l5| , Eq. (4) , where a static source expression 
for a black-hole binary has been derived for the next-to-leading order spin-squared terms: 



-^matter 

S^, static 2/7Z 



Am 



— (^61) + (7^1) + J-g mn jP^^ 7 n *,Sitf Wi 

mi V / mi v ' M 8mi F v 

L ( 7 ^ 7 m « 7 « m 5H„Si ifc 5i) . 



(4.74) 



with, in the present approximation, 1^ = m\^ k Q k ^ and Suj — tkijSik- Notice that contributions arising from the C3 
source term cancel each other, which is very nice, because C3 could not be determined and does also not contribute to 
the G 2 S 2 Hamiltonian. The explicit calculation then yields v\ = —2, v 2 = 3/4, v% = 3/4, 1/4 = 3/4, and v§ = —5/4. 
This result is fully consistent with the equations (|4.70p which were independently obtained by the Poincare algebra. 
Now all of the coefficients of the H S 2 p2 Hamiltonian and of the source terms have been fixed. They read 

m 2 3m 2 3 m 2 3m 2 , , 

at = - — 3-, a 2 = 0, a 3 = — 3, a 4 = --— 3, a 5 = 0, a 6 = -- — 3, (4-75) 
Am\ 8m^ 8 Am\ 

Ih = (3 3 = ft = ft = 0, ft = , ft = 7-^—, (4.76) 

4TOim 2 4mim 2 

3 9 3 3 3 15 

7i = T~2' 72 = --r— 2> 73 = 0, 74 = -— 3, 75 = _ ; 76 = -_ 77 = _^^ ) 4.77) 
4m^ 4m| 2mf mf 4m| 4mf 

^1 = 4! ^ 2 = 4' ^ 3 = ~24' ^ 4 = °' ^ 5 = 12' ^ 6 = 8' ^ 7 = 8' ^ 8 = _ 12' (4.78) 

It may be interesting to mention that in the ADMTT representation the A 4 source term is not present. Such a term 
emerges in [l6| with the values =F and with a 2 replaced by the square of the radius of the throat of a nonrotating 
black hole resulting from inversion symmetry at the throat. In Ref. [17] it was shown that a factor a 2 in the A4 term 
can also be generated by a specific deviation from the point-mass structure of a spherically symmetric body. 

V. THE TEST-PARTICLE LIMIT 

Next we will consider the test-particle limit of the energy of a black- hole binary up to the 4th order in l/ri 2 . To 
do this we plug into Eq. (5.2) of [18| the static ADM-gauged Kerr metric from Ref. Q labeled by particle number 1. 
The mass and momentum explicitly appearing in this equation will then get the (test-)particle number 2: 

/ \ ^l 2 

~Po = ~7i N u p 2j + Ni (m 2 + iiP2iP2j) ■ (5.1) 
The inverse 3-metric 7^ in ADMTT gauge up to the order 1/r 4 reads (G=l) 




ma 2 — 3m (a ■ n) 2 7m 2 a 2 9m 2 (a ■ n) 



2r 4 r 4 "V 'Hj 



Sij - h\( . (5.2) 



The result of this operation is 
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p\ TOlTO 2 3 TO1P2 1 ( 2 r TO l 2 



-Po =m 2 + ~ ~ h m 2 m 1 + 5 p 2 - 4roiai • (m 2 x p 2 

2m 2 ri2 2m 2 ri2 2rf 2 \ m 2 

1/1 o 25 TO? , 3 , .2 9 TOl ,/ \2 1 2 3 TOl 2 n 

m 2 m 1 - — p 2 + -TO1TO2 (ai • m. 2 ) +7 p 2 ( a i' n i2) ~ o m i m 2«i ~ 7 P 2 a i 



rj 2 \ 4 8 to 2 2 4to 2 2 4 to 2 

+ &m\a\ ■ (ni2 x p 2 ) ^ 

1/1 , 105 mf n 9 2/ \2 21 to? . 2/ \2 53 m? 9 2 

— o TO2m i + — p 2 - n^mi (ai • W12) + {p 2 ■ n 12 ) ai • n 12 ) - - p 2 (ai • ni 2 ) 

rj 2 \ 8 32 TO2 2 2 TO2 4 TO2 

7 to? . , 2 5 02 7 to? . , 2 2 23 to? 2 ? 21 o 

- 7 P 2 ' a i + ~m 2 m x a\ - - — (p 2 ■ n 12 ) a x + p 2 a 2 + — mjai ■ (n 12 x p 2 ) 

4 m 2 2 2 to 2 4 to 2 2 

+ 5mi (01 • n 12 ) 2 ai ■ (n 12 x p 2 ) - m\a^a\ ■ (n 12 x p 2 ) . 



(5.3) 



Now, from this limit we can also read off the coefficients f3 2 and P4 of Eq. (12.22[) in the ADM scheme, 

02 = -7^—. (5-4) 

4TOiTO2 

0a = 7-^—, (5.5) 

4TO1TO2 

in confirmation with the previously found ones. This is very strong evidence that we are on the right track with our 
method. 



VI. THE 2PN QUARTIC SPIN HAMILTONIAN 

To see that there is no S 4 Hamiltonian at the order 1/c 4 , we expand the Kerr metric in ADMTT coordinates up 
to the orders 1/r and a and show that such terms are not present at all and therefore cannot follow from a source 
with purely a terms. The strategy is the same as in [5|, but now we allow a terms and take the Kerr metric in 
quasi-isotropic coordinates of Eq. (43) in [5j] up to the aforementioned order and transform to ADMTT coordinates 
according to the formula 



with the 3-metric "ffj B ° 



4 DM = ir + ir^ 3 + ifk% + (6-i) 



qiso (s) . (2) (3) . (4)' la n 

Hi = l\i + Ta + Ta + Ta i I 6 - 2 ) 



iij ' ii] < H] • it] 



#' = # + ^^, (6-3) 



and the extended transformation vector 

laV \{a-n)a k 7 m 2 a 2 n k 7m 2 (a-n) 2 n k 



4 r 2 r 16 r 3 4 r 3 

1 a 4 n fc la 2 (a-n) 2 n k la 2 (a-n)a k l(a-n) 3 a k 
16 r 3 + 4 7 3 + 4 7 3 2 r 3 ' 



(6.4) 



We end up with a metric being the same as Eqs. (50) and (51) in [5j having all a 4 terms transformed away. This 
leads to the conclusion 

H S 4 = H s a = (6.5) 
for the quartic nonlinearities in the spin for the Hamiltonians linear in G. 
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VII. CONCLUSIONS 



The Hamiltonians we obtained are summarized as follows: 



G 

H s 3 lP i = -r 



at \ ( m2 C?2 5m 2 , c \ 

Si ' (ni2Xft) te Sr s (Si ' r,i2) 



(7.1) 



U SlS 2 pi -14" 



3 3 15 

2^2 S l s 2 ■ ("12 x p x ) + (Si ■ n 12 ) S 2 ■ (Si x p x ) - (5i • n i2 ) 2 S 2 • (m 2 x p x ) 



/ 3 



15 



n 12 • (Si x S 2 ) (^^2 Si ■ p-L - (Si ■ m 2 ) (Pi • TI12) 



(7.2) 



if 



G 
' 12 



4^3(^1 -Si) +^(P l '"i 2 ) Sj 



2 , ^2 2 ''""J 2 /o .. >2 



8m x 



8m x 



p : (Si ■ n 12 ) 



■"" ,2 (ft • n ia ) (Si • n 13 ) (p a • Si) - — L_ pls\ 



Ami 



\m\m 2 



Am\m 2 



3 9 3 

+ Irnf ' P2 ) S 'i - ^2 (Pi 'Pa) (Si -n 12 f - ^ (Pi ' "12) (p 2 • Si) (Si • n 12 ) (7.3) 



3 

H 2 (P2 ' "12) (Pi ' Si) (Si • n 12 ) 
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4m| 



4m^ 



(Pi • "12) (p 2 • "12) Si 



(Pi • "12) (p 2 • "12) (Si • ni 2 ) z 



plus the ones with interchanged indices (1 <-> 2), and 



H S 4 = H S 4 = . 

The corresponding sources in the constraint equations take the form 



(7.4) 



E 



b L 



1 

4m fc 



- ^<>."j - \p b • (a 6 x a) + ^Wy + ™ 2 ^ 1/2 - -^ n l 

1 „ „,o 1 



+ 7fr(Pb ' d )QbPbtd 3 



24m b 



(7.5) 



<7';' ( - j^PbkdA + ^p bt d k di + i (p 6 • a) <5 /4 <9 fc 

m 6 , „, 1 „,,., „ „ \ 1 



(7.6) 



The only missing Hamiltonians at the formal order of 1/c 4 are the G 2 Hamiltonians H S 2 + H S 2 resulting from 
gravitational nonlinearities of the Einstein field equations. These Hamiltonians will be presented in a forthcoming 
paper [151 ]. 
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